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The Kosterlitz-Thouless transition is studied from the representation of the systems's ground state wave func- 
tions in terms of Matrix Product States for a quantum system on an infinite-size lattice in one spatial dimension. 
It is found that, in the critical regime for a one-dimensional quantum lattice system with continuous symmetry, 
the newly-developed infinite Matrix Product State algorithm automatically leads to infinite degenerate ground 
states, due to the finiteness of the truncation dimension. This results in pseudo continuous symmetry sponta- 
neous breakdown, which allows to introduce a pseudo-order parameter that must be scaled down to zero, in 
order to be consistent with the Mermin-Wegner theorem. We also show that the ground state fidelity per lattice 
site exhibits a catastrophe point, thus resolving a controversy regarding whether or not the ground state fidelity 
is able to detect the Kosterlitz-Thouless transition. 
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Introduction. The Kosterlitz-Thouless (KT) transition (11], 
an infinite order transition from a critical phase to a gapful 
phase, is ubiquitous in quantum systems on an infinite-size 
lattice in one spatial dimension. It describes one of the in- 
stabilities of the Luttinger liquid induced by an marginal per- 
turbation. Remarkably, it does not arise from the long-range 
ordering, thus it falls outside of the conventional Landau- 
Ginzburg-Wilson paradigm, in which the most fundamental 
notion is spontaneous symmetry breaking (SSB) , with 
the symmetry-broken phase characterized by a nonzero lo- 
cal order parameter. Indeed, the type of internal order that 
a system possesses is profoundly affected by its dimension- 
ality. In quantum field theory, the Mermin-Wegner theorem 
states that continuous symmetries cannot be spontaneously 
broken in systems with local interactions in one spatial dimen- 
sion [4]. That is, in quantum systems with a continuous sym- 
metry in one spatial dimension, a genuine long-range order 
is destroyed by quantum fluctuations. Instead, a quasi-long- 
range order occurs in the critical phase, which is characterized 
by a power-law decay in correlation functions. 

The peculiarity of the KT transition makes it very difficult 
to map out the ground state phase diagrams for various quan- 
tum lattice systems in one spatial dimension. Indeed, a nu- 
merical analysis of the KT transition suffers from patholog- 
ical problems. One of these problems is that the finite size 
scaling technique jsj], which is successful for second order 
quantum phase transitions (QPTs) , can not be applied to 
the KT transition |8], since there are logarithmic corrections 
from the marginal perturbation. This motivated Nomura and 
Okamoto |@] to develop the so-called level spectroscopy to 
overcome these difficulties by combining the renormalization 
group, conformal field theory and the symmetry. This raises 
an intriguing question as to whether or not there is a unifying 
practical approach to different types of QPTs. 

In this Letter, we address this issue from a novel perspec- 
tive. The new inputs come from recent advances in both clas- 
sical simulations of quantum lattice systems and our under- 
standing of QPTs. First, a tensor network (TN) representa- 



tion of quantum many-body wave functions provides an ef- 
ficient way to classically simulate quantum many-body sys- 
tems, which include Matrix Product States (MPS) in one spa- 
tial dimension and the Projected Entangled-Pair States (PEPS) 
in two or more spatial dimensions. Second, two novel ap- 
proaches to study QPTs have be en prop osed in terms of both 
entanglement III and fidelity El El El H El El ■ In 
Refs. [12, ll], it was argued that the ground state fidelity per 
lattice site, combining with a practical way to compute it us- 
ing the TN algorithms for translationally invariant infinite-size 
systems, is able to capture the many-body physics underlying 
various quantum lattice systems in condensed matter. 

Although an MPS is well suited for the description of a 
gapped quantum state, it remains mysterious how it works for 
quantum states in the gapless regime. This casts doubts on 
the applicability of the NT algorithms to the study of the KT 
transition. Our approach is based on the observation that, in 
the gapless regime, the newly-developed infinite MPS (IMPS) 
algorithm automatically leads to infinite degenerate ground 
states, due to the finiteness of the truncation dimension. This 
results in pseudo continuous SSB iEtIi . which allows to intro- 
duce a pseudo-order parameter that must be scaled down to 
zero, in order to be consistent with the Mermin-Wegner theo- 
rem. We also show that the ground state fidelity per lattice site 
exhibits a catastrophe point 1 18] at a pseudo critical point for 
any finite truncation dimension. Normally, an extrapolation 
to infinite truncation dimension should be performed to de- 
termine the KT transition point. This resolves a controversy 
regarding whether or not theground state fidelity is able to 
detect the KT transition HslEai . 

Matrix Product States on an infinite-size lattice in one spa- 
tial dimension and continuous symmetries. The IMPS algo- 
rithm is a variational algorithm to compute the MPS repre- 
sentations of ground state wave functions for translationally 
invariant quantum s^^tems on an infinite-size lattice in one 
spatial dimension 11911 . Assume that the Hamiltonian takes 
the form: H - 2,- /i''-'^'^, with /jt'-'+'l being the nearest- 
neighbor two-body Hamiltonian density. The problem for 
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finding the system's ground state wave functions amounts 
to computing the imaginary time evolution for a given ini- 
tial state mO)}\: |*P(t)> = exp(-//T)|»I'(0))/| exp(-i/T)|»I'(0)>|. 
An efficient way to do this is based on the observation that 
the imaginary time evolution operator is reduced to a prod- 
uct of two-site evolution operators acting on sites / and i + I: 
U(i,i -t- 1) = exp(-/!t'''+'*J(5T), 6t « 1, as follows from the 
Suzuki-Trotter decomposition [20], and the fact that any wave 
function admits an MPS representation in a canonical form: 
attached to each site is a three-indices tensor or and 
to each bond a diagonal matrix Aa or Ag, depending on the 
evenness and oddness of the i-th site and the i-th bond, respec- 
tively. Here, s is a physical index, s = I, ■ ■ ■ ,d, with d being 
the dimension of the local Hilbert space, and I and r denote 
the bond indices, l,r - 1, ■ ■ ■ ,x, with;^^ being the truncation 
dimension. The effect of a two-site gate U{i,i + 1) may be 
absorbed by performing a singular value decomposition. This 
allows to update the tensors involved in the MPS representa- 
tion. Repeating this procedure until the ground state energy 
converges, one may generate the system's ground state wave 
functions in the MPS representation. 

It is known that the iMPS algorithm yields the best approx- 
imation to the ground state wave function (for a fixed trun- 
cation dimension x), as long as the ground state is gapful. 
However, it is surprising to see that it also works for quan- 
tum lattice systems with continuous symmetries in a gapless 
regime. The key point here is that continuous symmetries 
and the translational invariance under one site shifts cannot 
be maintained simultaneously during the imaginary time evo- 
lution 02 ill . Indeed, in order to mimic the gaplessness of ex- 
citations in the gapless regime, the iMPS algorithm resorts to 
the Goldstone mode. That is, the best approximation to the 
ground state wave function is not unique, if the translational 
invariance under one site shifts is maintained i22ll . Instead, 
infinite degenerate ground states are generated, each of which 
breaks the continuous symmetry. From now on, we restrict 
ourselves to the symmetry group t/(l). 

Pseudo continuous symmetry spontaneous breakdown and 
pseudo-order parameters. For a U{1) invariant quantum 
system on an infinite-size lattice in one spatial dimension, 
the iMPS algorithm automatically produces infinite degener- 
ate ground states, each of which breaks the t/(l) symmetry. 
Moreover, the symmetry breakdown results from the fact that 
an initial state has been chosen randomly. That is, aphe- 
nomenon occurs which shares all the features of an SSB 12 13]. 
Indeed, the implication of an SSB is two-fold: first, a system 
has stable and degenerate ground states, each of which breaks 
the symmetry of the system; second, the symmetry break- 
down results from random perturbations. In addition, such 
a "symmetry-breaking order" may be quantified by introduc- 
ing a local "order parameter", which may be read off from the 
reduced density matrix on a local area I23ll . However, this 
is in apparent contradiction with the Mermin-Wegner theo- 
rem, which states that continuous symmetries cannot be spon- 
taneously broken for quantum systems in one spatial dimen- 
sion 1,4J. To resolve this contradiction, one has to require that 



the "order parameter" must be scaled down to zero, when the 
truncation dimension x tends to oo. In order to distinguish 
them from a genuine SSB and a local order parameter, we 
introduce the notions of a pseudo SSB and a pseudo-order pa- 
rameter In this scenario, the Goldstone mode survives as gap- 
less excitations in the critical phase, and the KT transition is 
a limiting case of the conventional symmetry-breaking order, 
when the truncation dimension goes to oo. 

Ground state fidelity per lattice site. Suppose that a f/(l) 
invariant quantum lattice system undergoes the KT transition, 
when a control parameter A crosses a critical point Ac. Let 
us see whether or not the ground state fidelity per lattice site 
is able to detect it. The ground state fidelity per lattice site, 
d{Ai,A2), is defined as the scaling parameter, which charac- 
terizes how fast the fidehty F(Ai,A2) = \Ci'(A2W(Ai)}\ be- 
tween two ground states [^(/li)) and l^'jA^ )) goes to zero 
when the thermodynamic limit is approached il2lll3[l . In fact, 
the ground state fidelity F(Ai,A2) asymptotically scales as 
F(Ai,A2) ~ d{Ai, A2)^, with L the number of sites in a finite- 
size lattice. Remarkably, d(Ai,A2) is well defined in the ther- 
modynamic limit, and satisfies the properties inherited from 
the fidelity F(Ai,A2): (i) normalization d(A,A) = 1; (ii) sym- 
metry c/(/li, /I2) = d(A2, Ai); and (iii) range < d(Ai,A2) < 1. 

In the 1 /(1) symmetric phase, the ground state is non- 
degenerate [^24], whereas in the t/(l) symmetry-broken phase, 
infinite degenerate ground states arise. This implies that, if 
we choose ^(/la) as a reference state, with A2 in the U(l) 
symmetric phase, then the ground state fidelity per lattice 
site, d{Ai,A2), cannot distinguish infinite degenerate ground 
states in the U(l) symmetry-broken phase. This follows 
from the fact that <»I'(^2)l%(^i)) = mWim.i^i)) = 
(i'iA2)\^,j+^iAi)), for any large but finite size L, with U(^) 
being an element of the symmetry group U(l), and the sub- 
script 77 in \^',j{A)} labeling the eigenstates of the U(l) gener- 
ator. However, if we choose *I'(/l2) as a reference state, with 
A2 in the U{1) symmetry-broken phase, then d(A[,A2) is able 
to distinguish infinite degenerate ground states. Therefore, the 
phase transition point A^ manifests itself as a catastrophe point 
for any finite x- The catastrophe point disappears as x tends 
to 00. However, the critical point Ac may be determined by 
performing a scaling analysis with reasonably small ;^f's [i25J. 

Examples. The first example we shall investigate is the spin 
1/2 XXZ model described by the Hamiltonian: 



(1) 



where S „ {a - x, y, z) are the Pauli spin operators at the site 
/, and A denotes the anisotropy in the internal spin space. The 
model is in the critical regime (CR) for A e (-1, 1], with 
the same universality class as a free bosonic field theory, and 
exhibits anti-ferromagnetic (AF) and ferromagnetic (FM) or- 
ders, respectively, for A > 1 and A < - 1 . Besides a contin- 
uous t/(l) symmetry, the model also enjoys a Z2 symmetry, 
generated by the operation: 5 «-» S\} and — > -5'''. 
The second example is the spin 1 XXZD model with uniax- 
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Model 


a 




7 


6 


/J 


V 


xxz 

spin 1/2 





-0.00186 


0.17359 











0.16882 


0.04045 








XXZD 
spin 1 





-0.24659 


-0.27287 


0.00314 


0.08843 


-0.00204 





-0.24659 


-0.08133 


0.26049 


-0.07274 


-0.05033 



1.00 



TABLE I: The one-site reduced density matrices for the spin 1/2 

XXZ model (A = 0) and the spin 1 XXZD model with uniaxial 

single-ion anisotropy { J, = -0.5, D = -0.3): p^xxz = + °''^c'' 
1 , 



+ 

3 (1 - 2^6)7 + aS^ + piS^f + ^ySf + 

y!26Sf +n{{Sff - {Sff) + viSfsf + Sfsf). The two degenerate 
ground states originated from two randomly chosen initial states are 



connected via 51- 



['■]' 



: COS 6*5 sin (95 l'',5l,''' 



-sin 05';' + cos 6*5,", 



with e = 77.140 for the spin 1/2 XXZ model and 9 ■■ 
spin 1 XXZD model. 



72.001 for the 



ial single-ion anisotropy: 

H=Yj (5['l5f' i + 5M5f' i + y.5 W5f' i) + £ (5['l)2 



(2) 



where Sl^^ (a - x,y,z) are the spin 1 operators at the lat- 
tice site /, and J, and D are the Ising-like and single-ion 
anisotropies, respectively. For 7- = - 1 /2, there are two phase 
transitions, i.e., the KT transition from the large-Z) (LD) to a 
CR at D ~ 0.693, and a first order QPT from the CR to a FM 
phase at D 1.1 84, as D varies from oo to -oo [ 26] . 

Simulation results. In Table U we present the specific 
coeflicients in the one-site reduced density matrices for the 
spin 1 /2 XXZ model (A = 0) and the spin 1 XXZD model 
with uniaxial single-ion anisotropy ( 7- - -0.5, D - -0.3): 
Pxxz = 5^ + + ^^^'^ + y^y^ PxxzD = 5(1 - 2y6)/ -H 
a^I'l + j8(5")' + V2y5i'' + ^/26S^:^ + ^i(iS^^f - (^j/V) + 
v(5[''5y^ -I- Sy^S^^^). The two degenerate ground states origi- 
nated from two randomly chosen initial states are connected 
via = cosOSf + sm6Sf,Sf' = -sin^^l!'' + cosOSf, 
with 6 = 77.140 for the spin 1/2 XXZ model (the fidelity 
per lattice site between them is 0.9976) and 9 = 72.001 for 
the spin 1 XXZD model (the fidelity per lattice site between 
them is 0.9925). The IMPS simulations are performed for a 
randomly chosen initial state, with;i' to be 16. 

We may choose ((5^''), (5y')) as the pseudo-order param- 
eter for both the models [27]. In Fig. [r|(a) and (b), we plot 

the magnitude (9^ - -^JiS^^^}^ + (5y')^ of the pseudo-order 
parameter for the spin 1/2 XXZ model and the spin 1 XXZD 
model (with a fixed 7- - -1/2), respectively. The IMPS simu- 
lations are performed for a randomly chosen initial state, with 
the truncation dimension x to be 8, 16, 32, and 50. This 
indicates that a first order QPT and the KT transition occur 
at A = -1 and A = 1, respectively, for the spin 1/2 XXZ 
model, and that a first order QPT and the KT transition oc- 
cur at D = -1.184 and D = 0.827, respectively, for the spin 
1 XXZD model. Note that an extrapolation to x - has 
been performed for the pseudo transition points (at which the 
pseudo-order parameter becomes zero), yielding the KT tran- 
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FIG. 1: (color online) Main: (a) The magnitude of the pseudo- 
order parameter ((5,''), (5''')) for the spin 1/2 XXZ model and 
(b) that for the spin 1 XXZD model (with a fixed J- = -1/2). 
The IMPS simulations are performed for a randomly chosen ini- 
tial state, with the truncation dimension x to be 8, 16, 32, and 
50. Note that the pseudo-order parameter is discontinuous even 
at A = 1 for the XXZ model. However, this is an artifact 
of the IMPS algorithm. Inset: The pseudo-order parameter is 
scaled down to zero, according to = a^"*(l -t- ')> the 
critical phase, in the left-top inset in (a) for the XXZ model 
(with a = 0.7331, 0.8766, 0.9545, = 0.3947, 0.2327, 0.1296, c = 
0.7493,0.1898,0.0224 for A = 0.4,-0.4,-0.8, respectively), 
and the middle-bottom inset in (b) for the XXZD model (with 
a = 0.9454, 0.9725, 0.9887, i = 0.1212, 0.081 1, 0.0625, c = 
0.0405, 0.1314, 0.0322 for D = 0.4, -0.4, -0.8, respectively). More- 
over, an extrapolation to = 00 is performed for the pseudo transition 
points (at which the pseudo-order parameter becomes zero), yielding 
the KT transition point (D = 0.827) for the XXZD model in the left- 
top inset in (b). 



sition point {D - 0.827) for the XXZD model (see the left-top 
inset in Fig. [T](b)), which is larger compared to D ~ 0.693 
from the level spectroscopy ll26ll . The pseudo-order parame- 
ter is scaled down to zero in the critical phase, as shown in the 
left-top inset in Fig.[T](a) for the XXZ model, and the middle- 
bottom inset in Fig. [T](b) for the XXZD model, to keep con- 
sistent with the Mermin-Wegner theorem. 

In Fig.|2](a) and (b), we plot the ground state fidelity per 
lattice site, t/(Ai, A2) , for the spin 1/2 XXZ model li28] and 
the spin 1 XXZD model (with a fixed 7. = -1/2). We have 
chosen *I'(A2) {^{02)) as a reference state, with A2 {D2) in 
the U{\) symmetry-broken phase A2 = {D2 - -0.2), then 
t/(Ai, A2) (d(Di,D2) ) is able to distinguish infinite degener- 
ate ground states, with a (pseudo) phase transition point as a 
catastrophe point (which coincides with the pseudo transition 
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FIG. 2: (color online) The ground state fidelity per lattice site, 
d(Ai,A2), for the spin 1/2 XXZ model (a) and, d(Di,D2), for the 
spin 1 XXZD model (with a fixed /; = -1/2) (b). Here, the trun- 
cation dimension^ takes 8, 16, 32, and 50, respectively. We have 
chosen 'P(A2) ('P(Z)2)) as a reference state, with A2 (D2) in the 
(/(I) symmetry-broken phase A2 = {D2 - -0.2), then rf(Ai,A2) 
{d{Di,D2) ) is able to distinguish infinite degenerate ground states, 
with a (pseudo) phase transition point as a catastrophe point. Note 
that rf(Ai,A2) is discontinuous even at A2 = 1 for the XXZ 
model dl. 



point from the pseudo-order parameter for a given x, within 
the accuracy). Note that £/(Ai,A2) is discontinuous even at 
A2 = 1 for the XXZ model 12911 . However, this is just an ar- 
tifact of the IMPS algorithm, since all the degenerate ground 
states collapse into the genuine ground state as;^ — > 00. 

We have also computed the von Neumann entropy for the 
spin 1 /2 XXZ model and the spin 1 XXZD model (with a fixed 
J, = -1/2) (not shown here). It detects the QPTs [3^, but 
does not distinguish infinite degenerate ground states arising 
from a pseudo SSB for a givenx- 
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